Application of a minimum cost flow algorithm to the three-dimensional gauge glass 

model with screening 
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We study the three-dimensional gauge glass model in the limit of strong screening by using a 
minimum cost flow algorithm, enabling us to obtain exact ground states for systems of linear size 
L < 48. By calculating the domain-wall energy, we obtain the stiffness exponent 8 — —0.95 ± 0.03, 
indicating the absence of a finite temperature phase transition, and the thermal exponent i/ — 
1.05 ± 0.03. We discuss the sensitivity of the ground state with respect to smaU perturbations 
of the disorder and determine the overlap length, which is characterized by the chaos exponent 
( = 3.9 ± 0.2, implying strong chaos. 
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In high-temperature superconductors, which are 
strongly type-II, disorder plays an important role as a 
pinning mechanism for vortices. Without disorder, as 
has been suggested by Abrikosov within mean-field the- 
ory, flux lines form a triangular lattice. However, un- 
der the influence of an external current perpendicular to 
the field, the vortices will move since they experience a 
Lorentz force, and thus energy is dissipated destroying 
the superconducting state. It has been suggestedlilfi that 
disorder in the form of point defects may pin the vortices 
at random positions leading to a superconducting vor- 
tex glass phase, where the phase of the order parameter 
is random in space, but frozen in time, similar to spin 
glasses. 

A simplified model commonly used to study the vor- 
tex glass is the gauge glass, which is believed to be in 
the same universality class and contains the necessary 
prerequisites, disorder and frustration, for glassy behav- 
ior. Within this model, one neglects fluctuations in the 
amplitude of the order parameter and only considers the 
phase of the condensate. The Hamiltonian is given by 



H 



J ^ cos( 



(/.J - Ay - A/ fly) 



X a 



(1) 



where (pi is the phase of the order parameter on site i and 
J is the interaction strength, henceforth set to J = 1. 
The sum is over all pairs (i, j) of nearest neighbors on 
a simple cubic lattice of size N — . The quenched 
random vector potentials Aij are drawn uniformly from 
the interval [0, 27r] and represent the effect of disorder and 
an external magnetic field. Screening of the interactions 
between vortices is incorporated by the fluctuating vector 
potentials which are integrated over from — cxd to oo 
under the gauge constraint V • a = 0, and Ag denotes the 
screening length. The limit Aq — > corresponds to strong 
screening, whereas Aq ^ oo is the limiting case without 



screening. The last term describes the magnetic energy, 
and is the sum over all plaquettes of the lattice, where 
the curl is given as the directed sum of the around 
one plaquette. 

Most of the theoretical work so far has concentrated 
on establishing numerically the lower critical dimen- 
sion of the gauge glass model, both with and without 
screening of the interactions between vortices. With- 
out screening, there is no finite temperatiuje|-|transition 
to a vortex glass phase in two dimensionsuEl, whereas 
in three dimensions there is evidence for a finite Tc, as 
has been found bitirdSffii^iii wall renormalization group 



analyses (DWRG, 
Carlo simulations 



and finite temperature Monte 
, though due to hmited system siffg 
and insufficient statistics the earlier DWRG studiesoO'El 
could not fully rule out that the lower critical dimension 
is exactly d — ?>. 

Sufhciently close to the critical point, screening effects 
become important, since the correlation length ^ diverges 
more strongly than the screening length A -and the two 
length scales eventually become corruparablaJ. The effect 
of screening was investigaied in Ref.Q by a DWRG study 
and more recently in Ref.U by means of a finite temper- 
ature Monte Carlo simulation, and the results indicate 
that screening is a relevant perturbation, destroying the 
finite temperature transition in three dimensions, though 
the DWRG analysis could only be performed for rather 
small system sizes {L < A). 

In the present paper we reinvestigate the gauge glass 
model in the limit of strong screening by performing a 
DWRG analysis using exact ground states, which we ob- 
tain via a minimum cost flow algorithm from combinato- 
rial optimization. This algorithm allows us to study sys- 
tems with linear size up to L = 48, which is considerably 
larger than the system sizes in the previous studiesEltil. 
In addition, we study for the first time the sensitivity of 
the ground state configurations with respect to small pa- 
rameter changes, thereby obtaining the chaos exponent. 
We make use of the vortex representation of the gauge 
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glass model, which is obtained from the Hapultonian in 
Eq. (|l]) by making the Villain approximatior£3t3, which 
replaces the exponentiated cosine term in the partition 
function by a sum of periodic Gaussians, and then inte- 
grating out the spin wave degrees of freedom. Thereby 
one obtains the vortex HamiltonianH 



Hi 



b,)G(*-j)(J,-b,), 



(2) 



defined on the dual lattice, which again is a simple cu- 
bic lattice. The 3i are three-component integer variables 
running from — oo to oo living on the links of the dual 
lattice and satisfy the divergence constraint (V • J)i = 
on every site i. The bj are magnetic fields which are 
constructed from the quenched vector potentials Aij by 
a lattice curl, i.e. one obtains bj as l/(27r) times the 
directed sum of the vector potentials on the plaquette 
surrounding the link on the dual lattice bi lives on. By 
definition, the magnetic fields satisfy the divergence free 
condition (V • h)i = on every site, since they stem from 
a lattice curl. The vortex interaction is given by the lat- 
tice Green's function 



(27r) 



L3 E 



1 - exp(ik ■ (r^ - Tj)) 
2E„=i(l -cos(fc„)) + Ag- 



(3) 



In the case we are interested in, the strong screening 
limit Aq 0, G{i — j) reduces to G(0) = for i = j and 
G{i,j) = ji£27rAo)^ for i ^ j with exponentially small 
corrections^. Thus if we subtract J(27rAo)'^ from the 
interaction and measure the energy in units of J(27rAo)^, 
one obtains the simpler Hamiltonian 



Hv = 



We remark that Hy is not trivial due to the divergence 
condition (V • J)^ = 0. 

Finding the ground state of the Hamiltonian in Eq. 
subject to the constraint (V • J)i = can be inter- 
preted as a minimum cost flow problem in the language 
of combinatorial optimization. From this point of view, 
the problem can be restated as 



Minimize z(J) — ^^Ci(Ji) 



(5) 



subject to the constraint (V • J)i — 0, where the cost 
functions Ci(Ji) — (J^ — bi)^/2 have been defined. 

The algorit]jmi|-we shall use is the successive shortest 
path algorithnv^jiJl, which solves the problem in polyno- 
mial time in this specific case. For the implemeplation 
we made use of the LEDA programming libraryta. We 
were able to obtain exact ground states for system sizes 
up to L = 48 on ordinary workstations. The compu- 
tation time increases approximately with N'^, and one 
instance for L = 48 took about 2.5 hours computer time 
on a Sun Ultra2 (167MHz) workstation. We remark that 



a similar algorithm has already been used to calculate 
ground states of the solid-on-solid model for a surface on 
a disordered substrate, which also can be mapped on a 
minimum cost flow problcmED. 

Wc now present our results. First we address the ques- 
tion if the gauge glass model in the limit of strong screen- 
ing shows a finite temperature transition as presumably 
is the case without screening, i.e. we investigate the sta- 
bility of the ground state with respect to thermal fluc- 
tuations. We make use of the concept of domain wall 
renormalization, whieh has been applied to spin glasses 
in the same contextt2l. The idea is to study the energy 
AE necessary to flip a cluster on length scale L. For long 
length scales L, one expects that 



(6) 



where 9 is the stiffness exponent. The sign of 9 now 
determines whether there is a finite temperature phase 
transition. If 6 is positive, then the domain wall energy 
is increasing with cluster size, and one concludes that 
the ground state is stable with respect to thermal fluctu- 
ations. Correspondingly, if 9 is negative, the argument is 
that the ground state is unstable, since large cluster can 
be flipped by arbitrarily small energy. Thus in this case, 
the transition temperature is zero. 

The usual way to determine the defect- or domain 
wall energy is to measure the energy difference between 
ground states obtained for periodic and anti-periodic 
boundary conditions (be), respectively. For the model 
under consideration in the vortex representation (|^), 
however, one has to incorporate a_ boundary term to 
mimic the effect of a change in bojfi. This boundary 
term is rather inconvenient to implement in the mini- 
mum cost flow algorithm we are using, so we propose an 
alternative method to induce a domain wall or elemen- 
tary excitation. Note that the vortex Hamiltonian in Eq. 
(^) with periodic be and without the additional bound- 
ary term corresponds to fl-uetuating boundary conditions 
in the gauge glass modelcj'til. 

In the gauge glass one changes from periodic to anti- 
periodic be by adding tt/L to the vector potentials in one 
spatial direction, e.g. Afj Afj + ^. Such a shift of the 
vector potential, however, has no effect on the Hamil- 
tonian in the vortex representation (Q), since the mag- 
netic fields hi are constructed from the vector potential 
by a lattice curl and thus remain unchanged. Only the 
argument of the boundary term changes, which can be 
compensated by forming one or several vortex loops with 
total areafl'tl L^/2. From this observation we derive our 
basic idea: First we consider the vortex Hamiltonian Eq. 

with periodic be and calculate the exact ground state 
configuration {J*^}. The energy Eo{{3^}) of this state 
is obtained via Hy in Eq. (Q). We then determine the 
global flux / of this conflguration in one spatial direction, 
e.g. in the x-direction 



fx — ^ -^f ; 



(7) 
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where fx can be interpreted as a total winding number. 
Next, we gradually decrease all costs in the x-direction 
of the minimum flow problem in Eq. (||) simultaneously, 
which makes global flux in this direction energetically 
more favorable, whereas the cost for all topologically sim- 
ply connected loops (those with winding number zero 
around the 3d torus) remains unchanged. We reduce the 
costs until we obtain a configuration {J^} with global 
flux /a; + 1, which is an elementary low energy excitation 
with length scale L. We calculate the energy Ei of this 
configuration again with Hy in Eq. with the same 
magnetic fields b,; as before. The domain wall energy is 
then given by AE — Ei — Eq, which is always positive 
since the new state {J^} with flux fx + ^ corresponds to 
an excited state. 

In a small fraction (« 5%) of the samples, the flux 
changes discontinuously by more than one unit upon 
slowly decreasing the costs as described above. However, 
the resulting conflguration still represents an elementary 
excitation of length scale L and we can also use this con- 
figuration for calculating the energy. 

We note that it is easy to see that in the vortex 
(or Villain) representation of the pure classical three- 
dimensional XY- model which is given hy H = X^i^fj 
the elementary low energy excitation is indeed given by 
a configuration with one additional flux line of length L 
winding once around the 3d torus. The energy difference 
between the ground state {fx — 0) and the excited state 
{fx = 1) then simply is Ai? L, which corresponds to 
the energy of a spin wave excitation with minimum wave 
vector in the 3d XY-model in the phase representation. 



AE - with 



-0.95 ±0.03. 



(8) 
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FIG. 1. The domain wall energy [AiSjav in a log-log plot. 
The straight line is a fit to [A£;]av ~ with 6 = -0.95±0.03. 
This implies a thermal exponent of u = 1.05 ± 0.03. The 
disorder average is over 500 samples for L = 48, 1500 samples 
for L = 32, and for the smaller sizes several thousand samples 
have been used. 

Fig. 1^ shows the domain wall energy AE vs. L for 
L < 48 in a double logarithmic plot. One observes a 
straight line behavior which can be nicely fltted by 



Thus we reestablish that Tc = for the gauge glass 
model iBi-the strong screening limit, as has been found 
in Refs.BO. From the stiffness exponent 9 the thermal 
exponent which describes the divergence of the corre- 
lation length, can be calculated. For Tc = it is ^ ~ T~'^ 
and by equating the thermal energy with the energy of a 
low lying excitation on the length scale of the correlation 
length, it follows that 



1 



(9) 



From this relation we obtain v — 1.05±0.03, which agrees 
well with a result from a flnite temperature Monte Caiip 
simulation for the same model by Wengel and YoungO, 
who were able to study system sizes L < 12 and found a 
zero temperature phase transition with v = 1.05 ±0.1. 

Next we want to discuss the issue of chaos in the gauge 
glass model. From spin glasses it is known, that infinites- 
imal changes of parameters like the temperature or the 
couplings can have quite a dramatisreffect on the ground 
state or equilibrium configuratiorJ13c3. It is argued that 
a so called overlap length L* exists, which is a measure 
for the length scale, up to which the domain structure 
essentially remains unchanged after an infinitesimal per- 
turbation of a parameter. For length scales larger than 
the overlap length, the domain structure pchanges. For 
spin glasses, the overlap length is given bjii^ 



L* 



with C 



4 
2 



(10) 



where the result is obtained by equating the energy of a 
droplet excitation and the energy change caused by 
the perturbation which is proportional to L'^'l'^ , where 
ds is the fractal dimension of the droplets. 

To pursue a similar investigation in the gauge glass 
model we study the resulting changes in the ground state 
configuration when the random vector potentials Aij are 
perturbed by a small amount. To be specific, we define 
new vector potentials by A!^^ = Aij + , where is ran- 
domly drawn from the interval [—5, S\- We calculate the 
ground state for both realizations of the disorder {^y} 
and {^ij}, and define the distance D between the result- 
ing ground state configurations {J} and {J'} by 



(11) 



Note that for small values of Ds{L) the ground states are 
more correlated than for larger values. 

To determine the chaos exponent, we perform a scaling 
plot of the data for the overlap function Ds{L). Guided 
by Eq. (10), we attempt a scaling plot with the scaling 
variable LS^^'^, where C is the chaos exponent. Such a 
plot is shown in Fig. ^ and one observes that the data 
scales nicely with ( — 3.9 ± 0.2. This relatively large 
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value of C implies strong chaos, since in the limit of a 
vanishing perturbation (5 — > 0, the overlap length, which 
is proportional to 5^^^'^, increases slowly. The value of 
C ~ 3.9 ± 0.2 is considerably larger than for instance 
the one for the two-dimensional Isias spin glass, where 
^ — 0.95 ± 0.05 has been obtainedE3'E3. 




FIG. 2. Scaling plot of the overlap function Ds{L). The 
data scale nicely with the scaling variable L6^^^. The chaos 
exponent is given by ^ = 3.9 ± 0.2. The system sizes studied 
range from L = 4 to L = 32, where 5000 samples have been 
used for L < 20, 2000 samples for L = 24 and 500 samples 
for L = 32. 

Chaos cannot only be observed with respect to per- 
turbations of the disorder but also with respect to 
small temperature changesliB. For a continuous bond 
distribution— and low temperatures, the overlap length 
is expectedlia to behave as L*^ - T^^/C for = 0. 
The other relevant length scale is the correlation length 
^ ~ 5~'^, and for the 2d Ising spin glass with v Ki 2 and 
2/C w 2, the exponents characteriziag-the two length 
scales and ^ appear to be equalEJo, so one might 
speculate that they are related. However, for the gauge 
glass model we find 2/C ~ 0.5 and ly ^ 1, i.e. the diver- 
gence of the thermal correlation length is much stronger. 
It would be interesting to study chaos with respect to 
temperature changes in the gauge glass explicitly. 

Summarizing we studied the strongly screened gauge 
glass model in the vortex representation. Using a poly- 
nomial time minimum cost flow algorithm we could deal 
with much larger system sizes than considered before in 
the literature. We calculated the exact ground states and 
also configurations with a (global) low energy excitation. 
In this way we could perform a finite size scaling anal- 
ysis of the so called domain wall energy and obtained a 
pretty accurate estimate for the stiffness exponent that is 
6 = — 0.95±0.03. From this we can draw two conclusions: 
a) since it is clearly negative there is no superconduct- 
ing glass phase (or vortex glass phase) at non-vanishing 
temperature, and b) the thermal correlation length di- 
verges only with T ^ as ^th ~ T~'^ with an exponent 



v 1.05 ± 0.03, which is in agreement with Refs.Bllil. Fi- 
nally, we studied the effect of a small perturbation of the 
disorder on the ground state domain structure and found 
that the overlap length is characterized by the chaos ex- 
ponent ^ — 3.9±0.2. This is a pretty large value implying 
a fast destruction of ground state correlations by thermal 
fiuctuations. 
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